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r 

plenum 



hump speed 



yaw angle 

cavity beneath the wet deck of the XR-3 » 
formed by the sidewalls and the bow and 
stern seals, into which pressurized air 
is blown 

speed at which testcraft overcomes large 
low velocity drag and begins to operate on 
the cushion of cpatured pressurized air 
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I. INTRODUCTION 



A . BACKGROUND 

The Captured Air Bubble surface effect ship presents 
a new concept in shipping for the surface force of the 
U.S. Navy. Within a few years these vessels, supported 
partially by a cushion of air and partially by the more 
conventional hull buoyancy forces, will be joining the 
Fleet. In preparation for that event it is of extreme 
importance for detailed studies to be undertaken so that 
the capabilities and limitations of such vessels may be 
fully understood. 

The term surface effect ship applies to a broad spec- 
trum of different types of vessels. Among them are the 
Air Cushion Vehicle (ACV) , the Ground Effects Machine 
(GEM), the Semi-Submerged Ship (SSS) , the Captured Air 
Bubble (CAB), and others. This study is concerned with 
the CAB surface effect ship. 

The Captured Air Bubble craft consists of two rigid 
sidewalls which resemble catamaran hulls and flexible 
seals between the water surface and hull structure across 
the bow and stern of the craft. The enclosed volume is a 
plenum chamber into v/hich air is blown by a lift fan system. 
The pressure in the plenum chamber acting over the area 
of the plenum provides a lift force sufficient to support 
a large proportion of the total craft weight. The remainder 
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of the supporting force is made up by conventional hydro- 
dynamic forces acting on the sidewalls, plus some lift 
forces due to the bow and stern seals and aerodynamic 
lift at higher operating speeds. 

In the "on bubble" mode of operation, at the higher 
operating speeds, approximately two-thirds of the total 
craft weight is supported by the air cushion, and therefore 
relatively shallow sidewall immersion is required to support 
the remaining weight. The shallow draft produces a lower 
hydrodynamic drag and thus higher speed for a given power 
than could be obtained for a displacement hull. 

The XR-3 is a CAB surface effect ship of approximately 
three tons displacement. The XR-3 was built by the Naval 
Ship Research and Development Center (NSRDC) in 19 66. 
Following periods of testing by NSRDC, Annapolis Division 
NSRDC, and Aerojet-General Corporation, the XR-3 was de- 
livered to the Naval Postgraduate School (NFS) , Monterey, 
California in March 1970 for the investigation of basic 
and advanced surface effect ship technology. 

In addition to the XR-3 there exist other CAB test 
craft. Among these are the Aerojet-General 100-A and the 
Bell Aerospace Systems 100-B, both of approximately 100 
tons displacement. These two vessels were constructed as 
prototypes of a 2,000 ton ship. To date speeds in the 60- 
knot range have been attained by these vessels. 

Under a U.S. Government contract, Oceanics Incorp- 
orated developed a digital computer simulation of the 
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CAB ship loads and motions for the Surface Effect Ship 
Project Office (SESPO) in Washington D.C. This loads and 
motions program with the mathematical model of the Bell 
Aerospace Systems 100-B with six-degrees-of-freedom has 
been in use at the W.R. Church Computer Center at NPS 
since October 1972. 

In December 1973 » LCDR. D.G. Leo, USN and LT. R. Boncal, 
USN produced a mathematical model with six-degrees-of- 
freedom of the XR-3 test craft derived from the above- 
mentioned simulation model of the 100-B (Ref. 1). The 
modular construction of the 100-B simulation program 
facilitated the re-modeling of the program v'here necessary 
to fit the design differences between the 100-B and the XR-3- 

B. OBJECTIVES 

It is the objective of this thesis to validate the 
non-linear, six-degrees-of-freedom, mathematical model of 
the XR-3 CAB test craft in a comprehensive manner by direct 
comparison betv/een instrument-measured variables on the 
XR-3 and computed variables produced by the simulation 
programs. In the process of carrying out this objective, 
subroutine modifications were made to the simulation 
program to improve on the structure of the simulation model 
and justification is given in the sections that follow for 
the changes that were made. 

The data taken on the XR-3 measured in real-time, 
analog fashion using FM modulation was filtered, converted to 
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digital form, and processed in such a manner that it could 
be interfaced directly with the simulation program and a 
time history of true differences for each variable could 
be produced at each operating condition. 

The present study is limited to operation of the test 
craft in calm water conditions at the test site located 
at Lake San Antonio Reservoir, California. The range of 
operation of the test craft includes: 

1. Straight runs at various speeds with various 
locations of the center of gravity. The center 
of gravity was shifted by placement of two 125 
pound v/ater tanks at different locations on the 
deck of the test craft. 

2. Straight runs at various speeds in an unladen 
condition. 

3- Runs made at various speeds with various rudder 
angles applied to validate turning performance. 
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II. COLLECTION OF DATA 



A. INSTRUMENTATION ON THE XR-3 

In order to support the continuing studies being under 
taken at the Naval Postgraduate School on Captured Air 
Bubble vessels, the XR-3 has been instrumented with device 
capable of measuring and recording time histories of many 
different variables for later, post-run analysis (Ref. 2). 
The general configuration of the XR-3 is shown in Figure 1 
Table I shows a listing of the measurable parameters, the 
range of measurement, and the accuracy with which the 
measurements may be made. The locations of the measuring 
instruments are shown in Table II. The measured data are 
recorded on a Pemco Model 120-B magnetic tape recorder 
which features 14- tracks for recording data and an edge 
track for voice annotation during measurement runs. Each 
of the 14- tracks on the tape recorder is capable of record 
ing - 1 volt RMS with a tolerance of - 0.5$- 

B. POWER SUPPLIES 

The tape recorder is powered from a 110 volt 60 cycle 
auxiliary power unit, located on the aft weather deck as 
shown in Figure 1., through a Pemco a.c. to d.c. converter 
which provides a 26 volt d.c. pov/er source. The three 
accelerometers for measuring vertical, lateral, and fore- 
and-aft accelerations are powered from a 12 volt storage 
battery through a d.c. to d.c. converter which provides a 
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FIGURE 1. GENERAL CONFIGURATION OF XR-3 



Table I. XR-3 Measurable Parameters 



Parameter 


Simulation Program Range 


Accuracy 


Port Thrust 


Name 

THSTP 


0-500 lbs 


- 5 lb 


S tarb o ar d Thrus t 


THSTS 


0-500 lbs 


1 5 lb 


Velocity 


VEL 


0-40 kts 


- 1 kt 


Height 


None* 


- 2 feet 


- O.lin 


Bowseal Press. 


PBARB 


0-60 psf 


- 0.5 psf 


Plenum Press. 


PBAR 


0-60 psf 


- 0.5 psf 


Sternseal Press . 


PBARS 


0-60 psf 


- 0.5 psf 


Pitch Angle 


THETAR 


- 15 degs 


- 0.5 deg 


Roll Angle 


DPHI 


- 15 degs 


- 0.5 deg 


Yaw Angle 


DPS I 


- 180 degs 


- 0.5 deg 


Pitch Rate 


PDEG 


-30 degs/sec 


-0.5 deg/sec 


Roll Rate 


QDEG 


+30 degs/sec 


-0.5 aeg/sec 


Yaw Rate 


RDEG 


-30 degs/sec 


- 0.5 deg/sec 


Rudder Position 


DELRS 


- 45 degs 


not available 


Longitudinal 

Acceleration 


None 


- 0 . 2g units 


not available 


Lateral 

Acceleration 


ACCLAT 


- 0.2g units 


not available 


Vertical 

Acceleration 


ACCEL(3 ) 


- 0.8 to +1 . 2g 


not available 



^Height is used to compute Draft as explained in text. 



Table II. Location of Sensors 



Sensor 

Pressure 

Transducers 

Bow Seal- 

Plenum 

Stern Seal 

Velocity 

Height 

Gyros 

Pitch 

Roll 

Yaw 

Port Thrust 

Starboard Thrust 

Vertical 
Accelerome ter 

Lateral 
Accelerome ter 

Longitudinal 
Accelerome ter 



Distance Forward 
of Transom 


Distance from 
Center Line 


18'4" 






i 


12 ' 4| 


n 




S. 


1 ' 6 " 






i 


21'9" 

4’9i” 


down 

deck 


from 




24' 1" 

6" 


down 

deck 


from 





12 ' 6" 


2'4i” Port 


13' 


2'4|” Port 


12 ' 6" 


2’4i" Port 


0 


5' 6" Port 


0 


5' 6" Starboard 


13' 61" 


1 ' 9i" Port 


13 ' 6 " 


1'9" Port 


13' 5i" 


1 ' 8 - 2 -" Port 



19 



28 volt d.c. supply. The three gyros which measure pitch, 
roll, ya w, pitch rate, roll rate and yaw rate are powered 
by the 12 volt storage battery. The gyros also receive a 
1 volt d.c. excitation signal from the 110 volt a.c. power 
supply fed through an a.c. to d.c. converter. 

The height sensor mounted forward of the bow is 
powered directly from the 110 volt a.c. auxiliary power 
unit. The height measurement signal is converted to a 
-1.0 volt d.c. signal by a height sensor conditioner. The 
velocity probe, mounted on the same assembly as the height 
sensor, extends into the undisturbed water forward of the 
testcraft. The 28 volt d.c. power source for the velocity 
probe is the 12 volt battery through the d.c. to d.c. 
converter. The rudder position measurement is made by means 
of a load cell which receives a 12 volt d.c. supply from 
the battery. 

The two load cells for measuring port and starboard 
thrust and the three pressure transducers for measuring 
b ow seal, plenum, and stern seal pressures are all powered 
by a 28 volt d.c. supply from the 12 volt d.c. storage 
battery through the d.c. to d.c. converter. All five 
units also receive a 5 volt d.c. excitation signal from 
an amplifier- transducer package. 
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III. DATA PROCESSING 



A. ANALOG TO DIGITAL CONVERSION 

To provide measurements for validation of the mathe- 
matical model, data-gathering runs were made with the XR-3 
testcraft on four different days at the Lake San Antonio, 
California test site. The data was recorded on magnetic 
tape using FM modulation, then processed into a suitable 
form for comparison with corresponding values computed in 
the XR-3 simulation program on the IBM 360/67 Digital 
Computer. 

The data was converted from analog to digital form 
on an XDS-9300/CI-5000 hybrid computer at the NPS 
Computer Laboratory. From previous experience in data 
reduction by the XR-3 Research Group in the Aeronautical 
Engineering Department at NPS, it was found that filtering 
of the analog signal was required to remove unwanted high 
frequency noise from the recorded data. To provide this 
required filtering six Krohn-Hite digitally- tuned filters 
were assembled so that six tracks of recorded data could 
be processed simultaneously. 

B. DESCRIPTION 0? EQUIPMENT 

The analog tape containing raw data from XR-3 test 
runs was mounted on a Honeywell Tape Recorder/Reproducer, 
Model 96 for playback. Connections were made through 
trunk lines to lead the output from the six selected tracks 
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to the COMCOR Cl- 5000 Analog Computer as shown in Figure 2. 
Each channel of recorded data was then passed through an 
amplifier with a gain of 5 before being led out of the 
CI-5000 to one of the filters. The filters were operated 
in the lo w pass -maximally flat mode which provided a four- 
pole Butterworth response with an attenuation of 24db per 
octave from the cut off frequency. Cut off frequencies 
were set at 10 Hertz for all channels except port and 
starboard thrust which were set at 1.0 Hertz. Each 
signal at the output of its respective filter was led back 
into the CI-5000 for another amplification by a gain factor 
of ten. A Brush 8-Channel Recorder was also patched into 
the system at this point so that each of the six selected 
data tracks could be continuously monitored. Figure 2 
shows the signal flow patch through the hybrid computer 
system for analog-to-digital conversion. 

For analog-to-digital conversion, digital-to-analog 
conversion, and other applications, the CI-5000 Analog 
Computer and the XDS-9300 Digital Computer are intercon- 
nected through trunks. For A-to-D conversion in this 
application a sampling frequency of 200 Hz was chosen. In 
steady-state calm-water runs with the XR-3 none of the 
measured parameters can be expected to vary at frequencies 
greater than 10 Hz. Therefore, selection of a 200 Hz 
sampling frequency more than satisfies the requirements 
of the Sampling Theorem. 
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The first step in the A-to-D process was to zero all 
possible bias errors in the system. Calibration steps 
are recorded on separate channel strips on the analog 
data tape which indicate zero and maximum signal readings 
for each sensor. Bias potentiometers at the input of the 
first stage amplifier on the CI-5000 were adjusted to 
properly zero all channels during the pre-test run cali- 
bration of the XR-3 sensors. Each signal was monitored 
on a digital voltmeter during the calibration procedure. 

Execution of the FORTRAN program listed in Appendix A 
on the hybrid computer system enabled selective conversion 
from analog to digital form of specific time periods of 
analog data. Continuous monitoring of the voice edge track 
on the data tape provided a reference as to what the XR-3 
conditions were at all times. After each file of data had 
been digitized and written on the digital magnetic tape, 
the tape was re-wound and that file most recently recorded 
was run through the digital-tc-analog conversion. The 
D-to-A signal was monitored on the Brush Recorder and com- 
pared with the analog signal to verify the validity of the 
conversion process. The finished product of the A-to-D 
conversion was a reel of 7-track digital tape consisting 
of many files of digital data. Each file on the tape 
contained a time history of six selected variables measured 
for a specific operating condition on the XR-3- 
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c. INTERFACE WITH IBM SYSTEM 360/67 

Introduction of the 7-track tapes into IBM System 
360/67 at the W.R. Church Computer Center at NPS entailed 
some further processing of data. In the XDS 9300 a Binary 
number is represented by considering its octal representa- 
tion, made up of 3~t>it digits, whereas in the IBM 360 
the same number is represented by hexadecimal, made up 
of 4-bit digits . As an example of this difference , a 
decimal number (668) is (1234),-, octal and (29C)^£ 
hexadecimal. It was necessary, then, to write a FORTRAN 
program to read each record in a file from the 7-track 
tape, make the conversion from octal to hexadecimal, store 
the data in hexadecimal form on a temporary storage disk, 
and to print out a portion of each record in decimal form. 

A listing of the program is included as Appendix B. The 
program uses the pre-compiled assembler language subroutine 
FORM to convert the data to its hexadecimal representation. 
Figure 3 shows the flow of information through the IBM 
System 360/67 computer in this validation. 

The printed output from the conversion program is re- 
presentative of the data in its original measurement form 
multiplied by a gain factor of 50 in the A-to-D conversion 
process. At this point in the process the numbers printed 
out were compared with the signals recorded on the Brush 
Recorder during A-to-D conversion to insure that the 
accuracy of the data had been preserved. This was indeed 
the case with the exception of some small amplitude, 
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FIGURE 3, INFORMATION FLOW TMRUOGH IBM SYSTEM 360/67 



high-frequency fluctuations in the digitally reproduced 
data which did not appear to be present on the Brush 
Recorder tapes. 

The final phase in data processing was to recover the 
measured data into the units which correspond with the 
measurement, i.e. pounds of thrust, degrees of pitch angle, 
knots of velocity, etc. to facilitate direct comparison 
in the XR-3 simulation model. To accomplish the recovery 
of data it was necessary to read the converted data from 
the temporary storage disk, multiply each set of variables 
by the inverse of the appropriate scale factor used when 
recording the data on the testcraft, divide out the gain 
factor of 5° introduced in the A-to-D process, and store 
the result on a temporary disk for later access by the 
XR-3 simulation program. 

The height measurement made on the XR-3 during data 
collection runs does not have an immediately available 
counterpart in the XR-3 program as indicated in Table 1. 

A very close approximation to draft can be obtained from 
the height measurement through suitable scaling, however. 
The calculation of draft in the XR-3 simulation program 
is made v/ith reference to draft at the center of gravity 
of the craft. The height measurement is made at a sensor 
located 191" forward of the reference center of gravity 
and 6.375" down from the deck level of the test craft 
(see Figure l). The measurement of the height of the 
sensor above the surface of the water will vary as the 
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craft pitches up and down at the bow. Since pitch angles 
are, for the normal operating range, less than five degrees, 
the small angle approximation may be made for the sine of 
pitch angle. The draft measurement is synthesized by 
solution of the following relationship: 

DRAFT = SH + Z - HT + D*SIN Q 

where DRAFT = Water height above keel at the center of 
gravity 

SH = Height of sensor above CG = 7*1^5" 

Z = Distance of CG above keel = 30. 48" 

HT = Measured height of sensor above water surface 
D = Distance of height sensor forward of CG = 191" 
D*SIN0 = Variation of measured height due to pitch angle. 
A FORTRAN program was written to accomplish the process 
of re-scaling to measured unirs and to provide the close 
approximation to a draft measurement. A time index was also 
included in the program to facilitate plotting each variable 
versus time on a CALCOMP plotter. This feature proved to be 
valuable as a final step in the process because each variable 
could be closely examined in graphical form prior to acces- 
sing the data by the XR-3 simulation program. It was at 
this point in the procedure that some difficulties were 
encountered in early attempts at signal processing, as 
described in the following section. 

D. FREQUENCY COMPONENTS ENCOUNTERED IN DATA 

Examination of the CALCOMP plots of the digitally 
reproduced measurement data revealed significant small 
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amplitude, high frequency noise impressed on the measure- 
ment signal. The noisy signal was not immediately evident 
on the Brush Recorder tape, however on close inspection it 
was apparent that some noise was passing through the filters. 
Figures 4 , 5»6» and 7 are representative CALCOMP plots of 
processed measurement data of thrust, pitch angle, velocity 
and plenum pressure respectively. These graphs are in- 
cluded to show the effect of noise on the processed data. 

The source of noise present in the measurement data 
has not been established and can only be rationalized to 
be a combination of several possible environmental factors. 
First of these is vibration aboard the XR-3 testcraft from 
the twin engine propulsion system and fan drive system. 
Second, the testcraft is rarely run in completely glass- 
calm water conditions at the test site. Surface disturbances 
are nearly always present due to numerous pleasure craft, 
wind-induced ripples on the lake, and occasionally, the 
wake of the testcraft itself. Such disturbances can have 
a pronounced effect on the motion of the XR-3. Finally, 
60-cycle power supply interference can be assumed to be 
present due to the antenna effect of unshielded wiring 
both in the testcraft and in the NPS Computer Laboratory. 

1 . Elimination of Unwanted High Frequencies 

Recognizing that the noisy data rendered it all but 
unusable for any meaningful attempt at validation, methods 
were sought to eliminate the noise and to obtain a "clean" 
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signal and, at the same time, to preserve the validity of 
the actual measurement. The first effort was to filter 
at a lower frequency. The result was displeasing for, 
while the magnitude of the impressed noise was somewhat 
reduced, lower frequency oscillations which represented 
actual craft motions were nearly eliminated. The next 
possible solution was to employ two filters in cascade, 
thereby obtaining a low-pass filter with 48 db per octave 
attenuation above the cut off frequency. This procedure 
accomplished the main objective of removing virtually all 
the noise and passing only signals representative of the 
craft motions. The disadvantage of this method was that 
only three instead of six measurements were now possible 
because of the hardware limitations and synchronization 
became an insurmountable problem. 

The third possibility of removing the noise from 
the measurement signal was to proceed on the assumption 
that the noise was gaussian and white and to employ a 
smoothing or time -averaging technique. The technique em- 
ployed was to process the data as originally set forth and 
then, in the final scaling process on the digital computer, 
to obtain a moving time average. The algorithm implemented 
on the digital computer was to calculate the sum of the 
first ten points in the time history of a measured para- 
meter, divide the sum by ten and call the resultant figure 
the first point in the time history. Next, a counter was 
incremented by five and the process was repeated to obtain 
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the next point. In this manner the data is reduced from 
200 samples per second to 20 samples per second. The over- 
lapping smoothing technique tends to preserve low frequencies. 
By comparison there were insignificant differences (less 
than two per cent) between this smoothing process and the 
aforementioned cascade filtering. The advantage realized 
was a more economical use of available equipment (filters) 
and a great savings in time required for data processing. 

A listing of the program used for re-scaling the 
data to its original measurement form and including the 
smoothing routine is included as Appendix C. Figures 8, 

9, 10 and 11 are CALCOMP plots of smoothed measurement data 
of thrust, pitch angle, velocity and plenum pressure 
respectively. By comparison with the measurement data 
plotted in Figures 4,5»6, and 7> it is evident that the 
moving time average smoothing technique has preserved the 
measurement information of value and has eliminated the 
impressed noise. 
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IV. INTERFACE WITH XR-3 SIMULATION PROGRAM 



A. LOGIC ADDITIONS TO SIMULATION PROGRAM 

The XR-3 Loads and Motions simulation program was 
modified to the extent that, on execution of the program, 
the processed measurement data could be accessed if desired 
and on-line computation of differences could be made. The 
output from the program would then be a time history of 
the differences between computed and measured parameters. 
This entailed the addition of another option switch in the 
input data which is read by subroutine INCON. If the option 
switch is on or equal to one, subroutine INCON reads the 
measurement data into core from the disk where the data 
is stored. If the option switch is off or equal to zero, 
no reading of measurement data occurs . 

It was decided to restrict validation or data compari- 
son runs to ten seconds of simulation time because of core 
storage and computer time allocations. This decision set 
the size of the measurement data arrays at 200 elements. 

A COMMON block was added in subroutine INCON and RHS so 
that the measurement data could be available in subroutine 
RHS. As with other COMMON blocks, the feature to zero all 
elements of the added COMMON block was included in the 
main program. The following XR-3 test craft variables 
and the associated simulation names were selected for 
validation : 
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Total Thrust 



THST 



Velocity 


VEL 


Height 


DRFT 


Pitch Angle 


THE TAR 


Roll Angle 


DPHI 


Plenum Pressure 


PBAR 


Yaw Angle 


DEPSI 


Rudder Angle 


DELRS 


Vertical Acceleration 


ACCEL(3 ) 


Yaw Rate 


RDEG 



The necessary logic was added in subroutine RHS such 
that the appropriate element with respect to time in the 
measurement array could be compared with the computed 
variable and a difference formed. The resultant differ- 
ences were written into a temporary file for output by 
subroutine COLFIL after the simulation run reached the 
specified final time of ten seconds. Computation of 
differences was arranged such that positive values indi- 
cated that the computed parameter w as greater than the 
measured parameter for any given sample time. The print 
interval for output purposes v/as set at one point every 
0.05 seconds to correspond with the measurement data rate. 

B. METHOD OF THRUST INPUT AND OBSERVED DIFFERENCES 
The first step in the validation process was to 
initialize the mathematical model at a particular operating 
point with the expected trim conditions of velocity, pitch 
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angle, draft, and plenum pressure. The feature exists in 
the simulation program to provide a thrust map of up to 
25 points. Additional logic was included to select every 
fourth point in the 200 element array of measurement 
thrust to be applied to the mathematical model as an input 
variable and a ten second simulation proceeded to produce 
computed velocity. All other measurement parameters except 
velocity were ignored during initialization runs. The 
difference between computed velocity and measured velocity 
was taken as an output from the model. Following initial- 
ization of the model, another computer simulation was made 
using computed thrust as necessary to hold velocity con- 
stant. All measurement data for the run was read in and 
differences between computed and measured parameters were 
taken as output from the computer model. 

For validation studies using measurement data with 
very nearly constant thrust applied on the testcraft, the 
foregoing procedure could be shortened by including only 
the last phase. Provided that the initial trim conditions 
were known for input to the computer model, no loss of 
accuracy resulted in the process of validation. 

C. METHOD OF RUDDER AND THRUST INPUT TO THE MODEL 

Similar to the thrust map provision is a provision 
for a rudder map input to the computer model. Up to 25 
points may be included in a rudder map and the necessary 
logic was included in the program to select every fourth 



42 



point in the 200 element array of measured rudder angle 
for input to the mathematical model. The objective here 
was to apply the same inputs to the mathematical model as 
were used on the testcraft along with the initial trim 
conditions for each case. This method proved to be of 
little real value because a steering cable failure to the 
port propulsion engine on the testcraft rendered the rudder 
position measurement unreliable. To substitute for direct 
steering control on the port engine, the starboard engine 
was mechanically linked to the port engine by means of a 
long, rigid bar. Thus, the starboard engine controls the 
port engine in the rudder position and alignment of thrust 
vectors on the testcraft. Virtually all precision in rudder 
angle measurement is lost due to unavoidable play in the 
linkage and lack of direct steering control over both 
engines . 

Validation studies involving turning maneuvers with 
the XR-3 testcraft and the mathematical model were carried 
out by matching yaw angles and yaw rates for a spectrum of 
operating conditions and by noting the differences between 
computed and measured rudder angles for each operating 
condition. 
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V. PROGRAM SUBROUTINE MODIFICATIONS 



A. REPLACEMENT OF CENTER OF PRESSURE CURVE 

Newman and Poole reported in Ref. 3 the existence of 
a pressure wave and the attendant gradual deformation of 
the surface of the water due to a moving pressure distribu- 
tion with regard to CAB type hulls. NSRDC found that, in 
towing tank tests of CAB scale models, the craft tends to 
pivot in pitch motion about a point approximately 75~80fo 
of the wetted length of the sidewall forward of the tran- 
som (Ref. 4). Based on these findings, it has been 
postulated that the center of volume within the plenum 
chamber moves aft with increasing speed. 

Accordingly, the calculations made in the XR-3 sim- 
ulation program to account for a shifting center of pressure 
were removed from subroutines RKS and WAVES. Computations 
we re inserted in the program to account for the addition 
of an incremental wedge of volume to plenum volume which 
varies with craft speed. What follows is a description 
by subroutine of changes which were entered and the basis 
for the inclusion of the changes. 

1 . Subroutine INCON 

A previously unused space on the input data card 
number 00701 was used to read in the location of the pres- 
sure wave pivot point, measured in feet and tenths of 
feet forward of the transom. The statement XPWV = TEMP (3) 
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was inserted in Block 0? in the INCON subroutine. A value 
of 1?.2 ft. was used for XPWV in this study. 

The relationship for wave drag produced by the move- 
ment of the pressure bubble over the water surface is, 
from Ref. 5* 



''bubble = - (P. - P ,) C, 



4w 

/°gi 



The same formula is mechanized in the simulation program as: 



FXPWAV = -PBBAR * CF * PWVCON 
where CF = 0 .37 / ((U/FNCON) **1.5655981) and PWVCON = 

4.0* WEIGHT /(RHO *G* XLBW) . The coefficient 0.3? and the 
exponent 1.5655981 are correct for the length to beam ratio 
of the XR-3 as determined from Ref. 3* The variable name 
WATSLP is given to the slope of the water surface inside the 
plenum and is calculated as: 

WATSLP = PBBAR * CF * PWVCON / WEIGHT. 



WATSLP represents the angle at the apex of the wedge-shaped 
increment of volume to be added to plenum volume . Since 
the ratio of PBBAR*CF*PWVCON ( i.e. bubble drag) to WEIGHT 
will always result in a small number for speeds above hump 
speed, the small angle approximation may be used without 
loss of accuracy. Therefore, the volume increment due to 
the water slope to be added to plenum volume is given by 

0.5 * WATSLP * XL * AB 



where XL is the length of the plenum chamber and AB is the 
area of the plenum at the water surface. The quantity 0.5 
* WATSLP * XL is then the vertical cross sectional area 
of the wedge-shaped volume increment. 
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Two other calculations were added in subroutine 
INCON for use in other subroutines during the simulation. 

The distance from XPWV from the forward extremity of the 
plenum is calculated as 

XLXPWV = XLBW - XPWV 

where XLBW represents the plenum length at the water 
surface. Also the distance between the craft center of 
gravity and XPWV is calculated as 

XPWVXS = XPWV - XS 

where XS is the longitudinal center of gravity measured in 
feet forward of the transom. 

2 . Subroutine BOWSL 

In the computation for ELSKI , the wetted lenght of 
the bow seal, an incremental length is added to the wetted 
length of the seal. The increase is due to the water slope 
and is calculated simply as 

XLXPWV * WATSLP 

and added to the existing expression for ELSKI. 

3 . Subroutine STNSL 

Similar to the wetted length correction in BOWSL, 
an incremental correction is computed and applied in 
subroutine STNSL. Here the correction is opposite in 
sign and is applied to the wetted length of the stern seal, 
ELSKI as 

-XPWV * WATSLP. 

4 . Subroutine RHS 

The water slope correction is computed in subroutine 

RHS as a continuous variable over the period of any given 
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simulation time. Recognizing that WATSLP is calculated as 
the ratio of bubble wave drag to total craft weight and 
that bubble wave drag varies with velocity and plenum 
pressure, WATSLP must be considered as a continuous 
variable. Accordingly, the following computation was in- 
troduced in RHS : 

WATSLP = -FXPWAV / WEIGHT 

The plenum volume is corrected for the water slope on each 
pass through subroutine RHS by the addition of a volume 
increment as in subroutine INCONs 

+ .5 * WATSLP * XL * AB. 

5. Subroutine SIDEWL 

The addition of a wedge-shaped volume to plenum 
volume required an alteration to the computation of the 
drafts along the sidewalls inside the plenum chamber. 
Displacement of water is calculated as 

PBHEAD = PBAR / (RHO * G) 

where PBAR is the instantaneous plenum pressure. The im- 
mersion depth on the inside of each sidewall section is 
calculated in the relation 

DDIN = DD - WATSLP * (XPWVXS - XX ( ,K) ) 
where DD represents the uncorrected water level height 
above the keel along each sidev/all section and XX(J,K) 
is the distance from the transom of the center of each of 
the eleven stations along the sidev/all. 

A revision of the method of computing sidewall gaps 
at any given station was also necessary. A logical IF 
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statement was included to test each value of DDIN after 
it was computed. If DDIN is less than zero another logical 
IF statement must be tested to determine if a gap exists 
between the bottom of the sidewall (keel) and the water 
surface. If DSW ( J,K) , or the actual sidewall immersion 
depth, is less than PBHEAD , then a gap for that sidewall 
section is computed as 

GAP ( J,K) = - DDIN * ( 1 . - (DSW ( J ,K) )/ PBHEAD) 
where the computed GAP for any given station is in units 
of feet. 

Further background information of use in this 
study is contained in References 6-11. 

B. THRUST VECTORS IN SUBROUTINE PROP 

The propulsion system in the XR-3 consists of two 
forty horsepov/er outboard motors. The port screw has a 
right hand rotation and the starboard screw has a left 
hand rotation. Steering control is obtained by rotation of 
the motors in the horizontal plane. Thus, there is some 
rudder effect from the foil-shaped lower motor housings, 
however the predominant turning moment results from re- 
alignment of the propulsion thrust vectors when rudder angle 
is applied. In addition to the main thrust vector, directed 
in the line of craft motion for zero rudder angle, there 
exists a side thrust vector for each propeller. Since the 
propellers are counter-rotating, and due to the direction 
of rotation of each propeller, the side thrust vectors are 
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oppositely directed and, for zero rudder angle, are 
directed inboard. 

The assumption is made that both port and starboard 
propulsion motors rotate together through the same rudder 
angle. It is also assumed that the magnitude of the main 
and side thrust vectors are independent of rudder angle. 

The propeller thrust diagram shown in Figure 12 was 
used to verify and correct the force and moment equations 
for the XR-3 propulsion system. The figure displays the 
main thrust and side thrust vectors for both propellers 
positioned with right (positive) rudder angled, positive 
pitch angle©, and negative (port side down) roll angled. 
The thrust vectors act through the port and starboard pro- 
peller positions, XP, YP, and ZP from the craft center of 
gravity in the Body Reference coordinate system. 

The forces associated with the propeller alignment of 
Figure 1 are shown vectorially in Figure 13 in the Local 
Level coordinate system after having been corrected for 
pitch and roll angles. The assumption is made that pitch 
angle and roll angle are limited to small angles such that 
sin<9=©and cos© = 1, sin 0 = <P and cos ^ = 1. The follow- 
ing force equations result: 

T1X = THST1 * cos £ 

T2X = THST2 * cos S 
ST1X = - STHST1 * sin<£ 

ST2X = STHST2 * sinS 
T1Y = THST1 * sin S 
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FIGURE 13, VECTOR RESOLUTION OF PROPELLER THRUST FORCES 



T2Y = THST2 * sin S 
ST1Y = -STHST1 * cos 8 
ST2Y = STHST2 * cos i 
T1Z = THST1 * cos£ * & 

T2Z = THST2 * cosS 
ST1Z = STHST1 * cos £ * <t> 

ST2Z = STHST2 * cos S * 0 

Slimming the forces in each of the three corrdinate 
directions yields the following result: 

FXS = T1X - ST1X 
FXP = T2X + ST2X 
FYS = - T1Y - ST1Y 
FYP = T2Y + ST2Y 
FZS = - T1Z + ST1Z 
FZP = -T2Z - ST2Z 

The total force equations are found by summation of 
the force equations for the port and starboard propulsion 
systems in each of the three coordinate directions as 
follows : 

FX = FXP + FXS 
FY = FYP + FYS 
FZ = FZP + FZS 

Corrections were made to the force calculations of sub- 
routine PRO? from Reference 1 on the basis of the foregoing 
discussion . 
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VI . DISCUSSION AND EVALUATION OF RESULTS 



A. OBTAINING AGREEMENT AMONG VARIABLES 

Computer simulation runs in the initial stages of this 
study revealed that plenum pressure was consistently ap- 
proximately 2.5 psf higher than that recorded on the XR-3 
testcraft. Among other variables, thrust on the testcraft 
was higher than the computer model indicated for any given 
speed and calculated draft was greater by approximately 
three inches than the computer model produced. The decision 
was made to first obtain agreement in plenum pressure 
between the mathematical model and the testcraft. 

The bow seal leakage area input to the mathematical 
model was listed in Ref. 1 as one of the parameters not 
known to be exact. The estimate of leakage area was re- 
vised from 0.08 sq. ft. to 0.10 sq. ft. This slight alter- 
ation had the desired effect of reducing plenum pressure 
to agree with that measured on the XR-3 testcraft. An 
unexpected but welcome additional effect was nearly perfect 
agreement in thrust and draft. Aside from the alterations 
to the structure of the mathematical model there were no 
further changes made to the model during the course of 
the present study. 

1 . Straight Runs Made with Different C.G. Positions 

Validation runs were made interfacing data measured 
on the XR-3 testcraft with the computer model. The data 
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presented in Table III is a summary of average differences 
calculated between time histories of actual craft motions 
and simulated craft motions for the thrust, velocity, 
pitch angle, plenum pressure and draft variables. Also 
included is percentage differences for ease in qualitative 
comparison. It should be noted that percentage differences 
in pitch angle measurement appear to be quite large, how- 
ever in most cases the difference is within the measurement 
accuracy limitation of - 0.5 degree aboard the testcraft. 

The total craft weight on the measurement runs 
listed in Table III was 6270 lbs. In addition to the pilot 
there was a passenger in the co-pilot position and 250 lbs. 
of water in plastic containers aboard used for changing 
the location of the center of gravity of the XR- 3 - 

The data presented in Table IV is a summary of 
average time history differences with a fixed center of 
gravity on the XR- 3 - On the measurement runs listed in 
Table IV the total craft weight was 6050 lbs. Vertical 
acceleration was substituted for draft in this series of 
validation runs . 

Figures 1 4 , 15 > 16 , 17 . and 18 are typical of the 
graphical output from the XR -3 simulation program arranged 
to display the time histories of differences between com- 
puted and measured variables for this study. 

2 . Runs Made with Turns 

A series of eleven validation runs were made with 
various rudder angles at various speeds. A summary of the 
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TABLE III. SUMMARY OF DIFFERENCES IN STRAIGHT RUNS-A 
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TABLE IV. SUMMARY OF DIFFERENCES IN STRAIGHT RUNS-B 
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Craft total v/eight was 6040 lbs. in all runs. 



\ 




Figure 14. Plot of Thrust Difference Vs. Time 
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Figure 1 5 . Plot of Pitch Angle Difference Vs. Time 
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Figure 16. Plot of Velocity Difference Vs. Time 
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Figure 17. Plot of Plenum Pressure Difference Vs. Time 
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Figure 18. Plot of Draft Difference Vs. Time 
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average differences calculated bet ween actual craft motions 
and simulated craft motions during the turning maneuvers 
for thrust, velocity, pitch angle, roll angle, yaw angle, 
yaw rate and rudder position variables is presented in 
Table V. 

The relatively large differences in rudder angles 
in many runs is attributed to the lack of confidence in 
an accurate rudder position measurement. In the future 
when direct steering control of both propulsion motors is 
restored, it is recommended that confirmation be made of 
the results presented here with regard to rudder position 
and other variables in turning maneuvers. 

Figures 19 » 20, 21 and 22 are typical of the graph- 
ical output from the XR-3 simulation program which display 
the time histories of differences in roll angle , yaw angle , 
yaw rate, and rudder angle, respectively. 

B. LIMITED TO CALM WATER CASES 

It is recognized that much of the measurement data 
obtained at the XR-3 test site at Lake San Antonio, Cali- 
fornia contains small perturbations due to surface 
disturbances on the lake. No attempt has been made to 
measure the surface disturbances or their effect on the 
testcraft. Validation of the mathematical model in the 
presence of waves of any amplitude and wavelength is sug- 
gested as an area for further study. 



63 



TABLE V. SUMMARY OF DIFFERENCES IN TURNS 
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TABLE V. (Continued) 
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Figure 19. Plot of Roll Angle Difference Vs. Time 
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Figure 20. Plot of Yaw Angle Difference Vs. Time 
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Time 

(2 sec/inch) 



Figure 22. Plot of Rudder Angle Difference Vs. Time 
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A listing of the FORTRAN program for simulation of 
the XR-3 is included as Appendix D. A sample input data 
deck is also included. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 



This study has shown that the nonlinear mathematical 
model with six degrees of freedom adequately describes the 
motions of the XR-3 in calm water conditions. The variables 
used in the validation included thrust, longitudinal velo- 
city, pitch angle, plenum pressure, draft and vertical 
acceleration in straight runs. For validation of turning 
dynamics the variables included thrust, longitudinal velo- 
city, pitch angle, roll angle, yaw angle, yaw rate and 
rudder angle. Differences that were found to exist between 
the computer model and the testcraft were less than 5% for 
the majority of the variables compared and those differ- 
ences that exceeded 5^ are considered acceptable, given 
the range of the measured variables and the accuracy of 
the sensors aboard the testcraft. 

During the course of this study, it became apparent 
that the computer model exhibits a more oscillatory, or 
underdamped, response in pitch and roll motions than does 
the XR-3 testcraft. A concurrent study by LCDR. R.A. Finley, 
USN upon this aspect of the simulation with attention to 
the modelling of the bow seal and the lift fan system 
is reported in Reference 12. 

The method of data processing sot forth in this study 
to produce data in adequate form for interface with the 
mathematical model is recommended for use in future 
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validation efforts. A considerable amount of time and 
effort was spent in signal analysis and data processing 
to prepare the measured data in a form useful for compari- 
son with the mathematical model variables. In order to 
eliminate the requirement for smoothing the re-scaled 
data, it is recommended that a study be made using a lower 
sampling rate for the A-to-D conversion to see if the high 
frequency noise can be eliminated from the digitally re- 
produced data. It is also recommended that a vibration 
study be conducted aboard the XR-3 in an attempt to deter- 
mine if the source or sources of the high frequency noise 
is the result of internal ship vibration acting on the 
measurement transducers, gyros, and accelerometers. The 
current data acquisition system installed in the testcraft 
should be adequate for such a study. 

A very fertile area for future study would be valida- 
tion of the computer model in the presence of waves. 
Instrumentation exists to measure wave conditions at the 
XR-3 test site during measurement runs and the simulated 
waves could be input to the mathematical model. Direct 
comparison of craft behavior and simulation model calcula- 
tions could then be obtained to produce an evaluation of 
the accuracy of the simulation in sea states. 
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APPENDIX A 



LISTING OF FORTRAN PROGRAM FOR ANALOG TO 
DIGITAL CONVERSION OF MEASURMENT DATA 



DIMENSION I BUF( 40 96,2) ,LOCB(-l ’ l ) ,BADREC( 100) 

INTEGER RECNUM, BADREC 

NAMELIST NREC, NS AMP , NCHAN, I T APE , NOEL , NT 

1 I N PUT (101) 

NWGRDS=NSAMP* NCHAN 
LOCB( - 1 ) =L OCF ( IBUF( 1,1)) 

IF NREC.GT .1 ,LOCB( 1 ) =LOCF ( IBUF( 1,2)) 

IF NREC.EQ.1»L0C3( 1 ) = L OCB ( - 1 1 +N WORDS 
I ND = 2 
GO TO 15 

2 N B= 1 

R ECNUM = 0 
NEWBUF=LCCB(1) 

I F ( IABS ( I NO) .NE.2 ) GO TO 94 
CALL Wk ITECLOCK(O) 

CALL ADSTART (NCHAN, L0CB(-1 ) , NEWBUF , NS AMP , RECNUM, 1 1 S ) 
CALL MTRDYi I TA PE , LOCB ( - 1 ) ,L0C3( 1 ) ,NWORDS» IND) 

4 IF IND.EQ.l, GO TO 4 
GO T0( 4,3, 91 ,93 ) IND 

3 IF NT . E Q. 0 , GO TO 7 

CALL MTRDYK ITAPE + 1 , LOCB ( - 1 ) , LOCB( 1 ) , NWORDS , I NDJ 

6 IF IND.£Q.1,G0 TO 6 

GO T0( 6,7,91,93) IABS( IND) 

7 I F ( TEST 11) . GT . 0 ) GO TO 7 
I ND=2 

mo a r\— n 

I >1 'J M U — w 

CALL STARTCLOCK 
CALL ENAdLE 

5 CONTINUE 

10 GO TO 5 

11 N6=-NB 
NtWBUF=LOCB( NB ) 

IF( TEST ( 1 ) .GT.O. 0 R . RECNUM . G E. N REC ) CAL L D I SABLE 
GO T0(90, 12,91 ,92) IA3S( IND) 

12 CONTINUE 

I F ( IND. bT.O) CALL MTOJT(NB) 

I F( IND. LT.O >CAi_L MTOJT 1 ( NS ) 

IF(TEST(1) .LT.O. AND. RECNUM. LT.NREOGO TO 5 
CALL STOPCLOCK 
CALL ADSTOP 

CALL PROCESS ( I3UF , NSA MP , NCH AN+- 1 , 2 S ) 

OUTPUT ( 6 ) RECNUM 

IF NBAD-NE.O, WRITE! 6 ,106) ( B A DR EC ( I ) , I = 1 , N3 AD ) 

106 FORMAT ( $ BAD RECORDS ARE$,(/,I6)) 

15 OUTPUT ( 102) 'CPI ION=( I 1) ' 

READ ( 101 , 100 )NOPT 

100 FORMAT (ID 

GO TCI 1 ,2,30,40, 50,o0, 70 )NGPT 
30 OUTPUT ( 102 )’ EOF ON WHICH TAPE’ 

READ ( 101 , 100) N 
IF N. EO . 0 , GO TO 15 
ENDFI Lc (N) 

GO TO 15 

40 OUTPUT ( 102 )' REw IND WHICH TAPE’ 

READ( 101 , l 00 ) N 
IF N . EQ . 0 i GO TO 15 
R E W I N D ( N) 

GO TO 15 

50 OUTPUT( 102) ' SKIPF ILES=( 14) ' 

R E AD ( 101 , 101 )NF 

101 FORMAT! 14) 

OUTPUT! 102) ' ON WHICH TAPE' 
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READ( 101, 100)N 
IF N. EQ .0 , GO TO 15 
DO 55 1=1, NF 

51 CALL BUFFERING, 1 , IBUF< 1, 1) ,1, IN ) 

52 I F ( IN . L T . 2 ) GO TO 52 
I F ( IN . NE . 3 ) GO TO 51 

55 CONTINUE 

OUTPUT ( 102 )NF 
GO TO 15 

60 OUTPUT! 102) ' NUM WORDS TO LIST=(I4)« 

READ! 101 , 101 ) N W 

OUTPUT! 102) ’ FROM WHICH TAPE' 

READ! 101, 1 00 ) N 
IF N.EQ.0,GC TO 15 
601 IN =1 

CALL 3UFFERIN(.N , 1 , I BUF ( 1 , 1) , N WORD S , I N ) 

66 IF! IN . EQ. 1 )G0 TO 66 

62 GO T0(62»63»64,65) IN 

63 WRITE! 6, 102) 

102 FORMAT! 1H1) 

DO 631 1=1 ,NW,NCHAN 

WRITE (6 , 104) ! I BUF! J, 1 ) , J=I , H-NCHAN-1 ) 

104 FORMAT ( 12010) 

631 CONTINUE 

IF (SENSE SWITCH 1)601,15 

64 OUTPUT ( 102) ' EOF READ' 

GO TO 15 

65 OUTPUT ( 102) ' READ ERR’ 

GO TO 63 

70 OUTPUT! 102)'D/A FROM WHICH TAPE' 

READ! 101 , 100)N 

IF N . E Q . 0 » GO TO 15 

OUTPUT ( 102) • TYPE * C/R TO CONTINUE' 

INPUT! 101) 

77 IN =1 

CALL BUFFERING! ; 1 ? I BUF , NWORDS , 1 N ) 

76 IF! IN . EQ. 1 ) GO TO 76 

71 GO T0(7 1,72,6^,74) IN 

72 DO 73 1=1, NWCRDS 

73 IbUFII, 1)=IBUF(I,1)/2**10 
DO 75 1=1, NWORDS , NCH AN 

DO 750 J = 1 , NCHAN 
750 CALL DAC! J, IBUF! I+J-l , 1 ) ) 

M=NDEL 
CALL DELAY 
75 CUNT I NUE 

I F( SENSE SWITCH 1)77,15 

74 OUTPUT! 102) ' READ E RrsOR ' 

GO TO 7 2 

90 CALL DISABLE 
CALL AOSTOP 

OUTPUT! 102) ’RATE ERR ' , RECNUM 
GO TO 15 

92 IF NBAD.LT. 1 00,NBAD=NBAD+1 
BADREC ( NBAD) =RECNUM-1 

GO TO 12 

91 CALL DISABLE 
CALL ADSTOP 

OUTPUT! 102) - MT NOT READY - 
GO TO 15 

93 CALL DISABLE 
CALL AOSTOP 

OUTPUT! 102) - MT ERROR ON SPACING FROM LOAD POlNT- 
GO TO 15 

94 OUTPUT! 102) ’DELAY TIME BETWEEN RECORDS TOO SHORT’ 
GO TO 15 

END 
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SUBROUTINE PROCESS! I » N , NC » I R ) 
J = N*NC# 9 
CALL DELAY 

1 IF TEST ( 1 ) .LT. 0 » GO TO I 

IF TESTt2) .GT.O, RETURN IR 

RETURN 

END 
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SEVEN TRACK TAPE, CONVERT FROM OCTAL TO 
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